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Abstract 

In this paper, we extend the study of bremsstrahlung photon pro- 
duction in a quark-gluon plasma to the cases of soft static photons 
(q = o) and hard real photons. The general framework of this study 
is the effective perturbative expansion based on the resummation of 
hard thermal loops. Despite the fact that bremsstrahlung only comes 
at two loops, we find that in both cases it generates contributions of 
the same order of magnitude as those already calculated by several 
other groups at one loop. Furthermore, a new process contained in the 
two-loop diagrams dominate the emission of a very hard real photon. 
In all cases, the rate of real or virtual photon production in the plasma 
is appreciably increased compared to the one-loop predictions. 
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1 Introduction 



We consider the production of a real photon or of a lepton pair in a quark- 
gluon plasma. The plasma is assumed to be in equilibrium at temperature 
T. The theoretical framework used in the calculation is that of thermal 
field theory improved by the hard loop resummation [1 — 6] of Braaten and 
Pisarski: in this approach one distinguishes hard momenta, of order T, 
from soft momenta, of order gT, where g is the Quantum Chromodynamics 
(QCD) coupling constant assumed to be small (g <C 1). After resummation 
of hard thermal loops, one is led to an effective field theory from which 
observables can be evaluated perturbatively. 

The production rates of real or virtual photons have already been evalu- 
ated, at the one loop level, in the effective theory [7 — 11]. Concerning soft 
virtual photons, it was found that the rate of production is considerably 
modified and enhanced compared to the result of the bare theory. Besides 
the usual quark-antiquark annihilation process, there appear many produc- 
tion mechanisms, in particular processes where the photon is radiated off a 
(hard) quark in a scattering process where the quark is backward scattered 
in the plasma via soft quark exchange. This is to be contrasted to the result 
obtained in a semi-classical approximation [12 — 16] where the photon is ra- 
diated off fast quarks in scattering processes mediated by a gluon exchange: 
we call such processes bremsstrahlung emission of a photon. In this study we 
reconcile the two approaches and show that the bremsstrahlung processes 
favored by the semi-classical approximation appear at the two-loop level in 
the effective theory and that, in fact, they contribute at the same order in 
gT as the processes in the one-loop effective theory. Such a result should 
not be a surprise. 

Consider the case of a soft virtual photon. The rate of production is 
related to the imaginary part of the vacuum polarization diagram [17, 18]. 
In the one-loop approximation of the effective field theory, it involves only 
effective fermion propagators and effective vertices [7]. A dominant contri- 
bution to the rate arises when the internal fermion momentum is soft and, 
therefore, all effective propagators and effective vertices have the same order 
of magnitude as their bare counterparts. A close examination of the final re- 
sult shows, however, that it has a logarithmic sensitivity to scales of order T 
(see section 3.3.2): this means that such a diagram also receives a dominant 
contribution from hard fermion momenta. When the momentum becomes 
large, the hard thermal loop (HTL) corrections to propagators and vertices 
are suppressed by, at least, a factor g with respect to their bare counter- 
parts. This suppression factor g can easily be compensated by the larger 
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phase space available to a hard fermion (O(T)) compared to soft fermion 
(O(gT)), thereby leading to a contribution of the same order of magnitude 
from the soft region and the hard region of phase-space. Now, when an 
observable is sensitive to the thermal corrections of hard vertices and prop- 
agators, it is obvious that all such corrections should be taken into account 
for the calculation to be complete. Some of these thermal corrections are 
naturally included in the lowest order of the effective theory via the resum- 
mation of hard thermal loops. But, even if the HTL approximation is correct 
for soft external particles, it does not account for all thermal corrections to 
hard vertices and hard propagators. For instance, neglecting the external 
momenta as one does in the calculation of the hard thermal loops is no 
longer justified when these momenta are not soft. Besides, equally impor- 
tant may be the contribution arising from soft gluons in the loop giving the 
HTL when the external momenta are hard, due to the Bose enhancement 
of the soft gluon term. Within the effective theory, both types of additional 
thermal corrections to a hard propagator or vertex are taken into account 
by considering a one-loop correction to this propagator or vertex. 

In the calculation of the virtual photon production rate in the effective 
theory, soft gluon exchange appears in two-loop diagrams. It will be seen 
that the bremsstrahlung production mechanism is precisely given by these 
diagrams when the exchanged gluon is space-like. The evaluation of these 
diagrams is discussed below. These contributions are clearly not included 
in the effective one-loop diagram. This is obvious when looking at the cor- 
responding physical processes and it manifests itself in the result by the 
calculated rate being proportional to the square of the thermal gluon mass 
nig, in contrast to the one- loop result where only the thermal quark mass 
appears [7]. Another important contribution of the two-loop diagrams is 
associated with time-like gluon exchanges: physically this represents QCD 
Compton scattering and quark-antiquark annihilation to produce a gluon 
and a photon. To evaluate this properly requires care since the two-loop 
diagrams with a hard time-like gluon exchanged are already part of the one- 
loop diagram with effective propagators and vertices. Taking into account 
the contribution of the appropriate counterterms in the effective lagrangian 
will prevent double-counting and allow the correct evaluation of the soft, 
time-like, gluon contribution. 

The case of soft real photon production follows essentially the same pat- 
tern, except for the crucial fact that, the external line in the vacuum po- 
larization diagram being massless, collinear divergences appear when eval- 
uating the two-loop diagrams: the quasi-overlap of two such divergences, 
associated with the fermion propagators, leads to an enhancement factor 
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of 0(1/ g 2 ) [19,20]. The paradoxical result then follows that the one-loop 
contribution is relatively suppressed by a factor g 2 compared to the two- 
loop one! The latter is entirely dominated by the bremsstrahlung process 
since the kinematical constraints require the gluon to be space like for the 
enhancement factor to occur. The Compton and annihilation mechanisms 
are sub-dominant. 

The case of a hard real photon, of momentum of 0(T), shares features 
with both cases above. The one-loop approximation has a logarithmic sen- 
sitivity to the hard fermion momentum in the loop [10, 11]. The two-loop 
bremsstrahlung has a I/5 2 collinear enhancement, as in the real soft photon 
case, which however is compensated by a factor m 2 /q 2 ~ g 2 , where q is the 
photon momentum, leaving the bremsstrahlung contribution at the same 
order in gT as the one-loop contribution. 

In the following we are concerned mainly with the bremsstrahlung part 
of the two-loop diagrams, and leave the discussion of the Compton and anni- 
hilation processes and their interplay with the counterterms of the effective 
theory to future work. We do not discuss the production of soft real photons 
since this has already been studied in detail in [19]. 

In the next section we derive the general expression for the (real or vir- 
tual) photon production rate at the two loop level. Then we consider the case 
of soft virtual photons produced at rest in the plasma and derive the leading 
behavior analytically. We compare to the one-loop results and show that the 
bremsstrahlung contribution is numerically dominant although both contri- 
butions are technically of the same order in g. The semi-classical approach 
is then discussed and it is shown that even though the approximations inher- 
ent in the semi-classical approach are not really justified in the case of soft 
photon production in a quark-gluon plasma, it leads to a result quite com- 
parable to that obtained in thermal field theory. Turning to the case of hard 
real photons, it is shown that the bremsstrahlung mechanism is of the same 
order as the already calculated one- loop result. Carrying out a more detailed 
comparison with the latter approach, it is found that the bremsstrahlung 
process dominates over the one-loop result for photon momentum of 0(T) 
but is relatively suppressed by a logarithmic factor for hard enough pho- 
tons. For very hard energies, we find that the photon production is in fact 
dominated by a new process consisting of a qq annihilation where the quark 
or antiquark undergoes a scattering in the medium. We summarize all the 
thermal field theory results concerning real and virtual photon production 
in a concluding section. 

The role of counterterms in the application of the effective theory up to 
two- loops is discussed in a first appendix where the problem is also illustrated 
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in a simple example. In a second appendix, the importance of phase space 
factors in thermal calculations is emphasized and the difference with the 
zero temperature phase space is made clear. 



2 Bremsstrahlung in thermal field theory 
2.1 Topologies involving bremsstrahlung 

Let us first recall the relationship between the photon production rate and 
the imaginary part of the retarded polarization tensor of the photon, as given 
by thermal field theory (we follow the notations of [21]). For real photons, 
this relation gives the number of photon emitted per unit time and per unit 
volume of the plasma as [17, 18]: 

dN do j , ra it/ \ 

2n B (g )Imn ^(q ,q), (1) 



dtdx (2vr) 3 2g D 

whereas for the production of a photon of invariant mass Q 2 decaying into 
a lepton pair we have: 

dN - dqodq "n.Wm^u). (2) 



dtdx 12vr 4 Q 2 

Basically, the above two formulae differ only by the allowed phase space for 
the photon, by an extra QED coupling constant when the photon decays into 
a lepton pair, and by the propagator of a heavy photon. It is worth recalling 
that these relations are valid only at first order in the QED coupling a since 
they do not take into account the possible re-interactions of the photon on 
its way out of the plasma nor the simultaneous emission of more than one 
photon. Nevertheless, they are true to all orders in the strong coupling 
constant a s . This should not be a serious limitation from a practical point 
of view since a <C a s . 

Let us now examine in which topologies the bremsstrahlung can ap- 
pear. It is worth recalling at this point that the retarded imaginary part in 
Eqs. (1) and (2) can be expressed as a sum over possible cuts through the 
diagram [22 — 24]. Therefore, we need to look for diagrams that will give 
bremsstrahlung processes once cut. A simple inspection of the processes in- 
volved in one-loop contributions (see [7, 19] for instance) shows that bremss- 
trahlung does not appear at this order. To see bremsstrahlung processes, 
one should consider the two- loop contributions of Fig. 1. The diagrams 
have been obtained via a strict application of the Feynman rules of the ef- 
fective theory [3,5], giving a priori effective vertices and propagators and 
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Figure 1: Two- loop contributions involving bremsstrahlung processes. A black dot 
denotes an effective propagator or vertex. Crosses are HTL counterterms. 

diagrams with counterterms 1 in order to avoid any double counting of ther- 
mal corrections already included at the one-loop level via the resummation 
of hard thermal loops, as outlined in appendix A. To make the connection 
with previous works [19, 20] easier, we mention that looking at two loop di- 
agrams in the effective theory is just a more rigorous way of doing what we 
might call "calculating one-loop diagrams beyond the HTL approximation" . 
Our present formulation is indeed more rigorous since it takes care of the 
counterterms, and also more positive since it does not assume a priori that 
one needs to go beyond the effective theory. Among all the possible cuts 
through the diagrams, those that correspond to bremsstrahlung necessarily 
cut the gluon propagator. Moreover, if L is the gluon 4-momentum, only 
the Landau damping part (L 2 < 0) gives bremsstrahlung, the L 2 > part 
rather giving Compton effect or q^-like annihilations [25]. There is another 
reason why the region I? > deserves a separate treatment: in the L 2 < 
kinematical domain, it is obvious that we cannot have contributions coming 
from the HTL counterterms since these counterterms involve only bare gluon 
propagators that don't have any imaginary part in the space-like region. On 
the contrary, in the L 2 > region, one should pay special attention to the 
counterterm diagrams in order to avoid any double-counting. Indeed, when 
the gluon becomes hard, we have a hard loop that may reproduce what 
is already included in the one-loop diagram via the effective vertices and 

1 These counterterms are nothing but the HTL contribution to the two or three-point 
function, with the opposite sign. Formally, they are necessary because one wants the 
effective theory to be just a reordering of the bare perturbative expansion, with the same 
overall Lagrangian. 
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propagators. From now on, we limit ourselves to the region where L 2 < 
and only to the true two-loop diagrams, leaving the 1? > region and the 
discussion of counterterms to future work. 

Moreover, since our main focus is on bremsstrahlung, we must retain 
from the cut quark propagators only the pole part and reject the Landau 
damping part, which would correspond to a different physical process. For 
the same reason, the cut should avoid going through an effective vertex 2 . 

2.2 General expression of two-loop contributions 

In order to obtain a bremsstrahlung contribution of the same order of magni- 
tude as the already calculated one-loop contributions, we need a hard phase 
space for the quark circulating in the loop, as explained in appendix B. We 
will see later that we have such contributions in the diagrams depicted in 
Fig. 1. Therefore, at leading order, we can use "bare" 3 vertices and propaga- 
tors everywhere except for the gluon propagator since the gluon can be soft. 
The diagrams we have to consider are therefore the simplified versions of the 
previous ones represented on Fig. 2. In the same figure, we have depicted 



Figure 2: Simplified two-loop contributions involving bremsstrahlung processes. The 
circled vertices correspond to the framework of the cutting rules. 

the relevant cuts as well as the arrangement of circlings that enables one to 
calculate the corresponding contribution in the framework of the thermal 
cutting rules. We have checked that the two cuts represented form a gauge 
independent set of terms, to which one should add the symmetric cut for the 
vertex diagram and a third diagram with the self-energy correction on the 

2 These extra requirements are not a claim that other configurations of the cut cannot 
give important contributions as well, but are dictated by our choice of looking only at 
bremsstrahlung. 

3 As in [19, 20, 26], we may have to keep an asymptotic thermal mass even in the hard 
region to regularize collinear divergences encountered when the external photon is on- 
shell. But contrary to [26], we don't need to add this mass by hand, since it is naturally 
contained in the hard limit of effective quark propagators which are our starting point in 
Fig. 1. 
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upper quark line. Since these two other terms give the same contribution 
as the previous two, we simply take them into account by multiplying the 
final result by an overall factor of 2. 

A straightforward application of the cutting rules valid for the "R/A" 
formalism, with the notations of [24], gives for the vertex correction: 



Imn ^{q ,q\ 



_Aii NC„ 

-imn /(?o,g)| wt « = -^ L 



d A P r d 4 L 



(2vr) 4 7 (2tt) 4 



ARR , „ „ ARR. , 

x e (Q, P + L,-R- L)g (-P - L, P, L) 

x g AAR (R + L, -L, -R)e AAR (R, -P, -Q) 

x Tr \rS-(P + L)YS A -(Ph, S RA (Rh°S- A (R + L)} D^(L) 



NC 



F 2 2 

eg 



d A p r d 4 L 



r d^P f 

J Wiv 



(27T) 4 / (2vr) 4 



KAL)-A(L) 



S (P)-S (P) [S (R + L)-S (R + L) 
x (n F (r ) - n F {p )) (n B (l ) + n F (r + l )) 
x S{R)S(P + L)P T pa L {L) Trace^| vcrtcx , 

where, following [19], we denote the fermion propagator: 



(3) 



S R ' A (P) = ?S R ' A (P) with P = (p , ^+M^p) (4-a) 
- =2^— = D2 - (4"b) 



P 2 ± i Po e P 2 -Ml± i Po e ' 
and the effective gluon propagator in a linear covariant gauge: 



D% (L) ee P PCT (L)A;' (L) + P pCT (L)A^ (L)+£L p L CT /L 2 (5-a 



L2-n TL (L) 

.2 



and p TL (L) = DiscA Ti (L) (5-b) 



n T (L) ee 3m; 



x" x(l — X ) , 

h — -In 

2 4 



n L (L) ee 3m 2 (l -x 2 ) 



I- -hi 
2 



X+ 1 

x- 1 

X + 1 

x — 1 



(5-c) 
(5-d) 



with P pa the usual transverse and longitudinal projectors in linear covari- 
ant gauges [2,27 - 30], M£, ee g 2 C F T 2 /4 [26] the asymptotic thermal mass 
of the quark, and m 2 = g 2 T 2 [N + N F /2]/9 the soft gluon thermal mass. In 
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this formula, e is the electric charge of the quark and therefore depends on 
its flavor. Likewise, we obtain for the second diagram: 

NC B 



Imn flA /(g D ,qr)| 



self 



-imn 



d A p 



d L 



x e (Q, P, - 

AAR I T-> 

x g (R + L 
x Tr [y<S A -(P) 7M <S 



(R,L, 



-R 



-L, 



„. AAR , „ 

i2)e (12, 



•P,- 



Q) 



F 2 2 

eg 



A (L)-A (L) 



(2vr) 4 7 (2tt) 



D P -{L) 



r d A P 


r d 4 L 


1 (2vr) 4 J 


(2vr) 4 



S R (P) 



S A (P) 



S (R + L)-S (R + L) P pa (L) Trace pcr 



x (5(P)) 2 (n F (r ) - n F (jp )) (n B (l Q ) + n F (r Q + l )) . (6) 

In the previous formulae, a factor S(R) without any R or A superscript 
simply denotes the principal part of the propagator. In other words, for these 
factors, the R or A prescription is irrelevant because the corresponding delta 
function is incompatible with the other delta functions present and therefore 
vanishes. 

We may notice the similarity between Eqs. (3) and (6). In particular, the 
same combination of statistical weights appear in both formulae, while the 
expressions in the square brackets simply express the cuts on internal lines. 
Moreover, when plugged into Eq. (1) in order to obtain the production rate, 
the sum of Eqs. (3) and (6) gives the more intuitive Eq. (62) in appendix B 
(to which one should add similar terms to take into account all the processes 
included in our diagrams). Therefore, the "R/A" formalism appears just as 
an efficient method to reorder the factors of the integrand in order to make 
it more compact and more convenient for the subsequent integrations. The 
drawback of this formalism is that it generates less intuitive expressions. 



2.3 Common part of the calculation 

The calculation of the Dirac's traces is of course common to both cases. We 
obtain for the self-energy insertion: 

Trace pCT Uf w -4 U~R 2 Q p R a - 4Q 2 R p R a 

-g pa (R 2 (R 2 -Q 2 ) + 2R 2 Q- L-2Q 2 R- l)1 , (7) 

and for the vertex correction: 

Traced | vertex w -4 \2R 2 P p Q a - 2{FTL) 2 R p Q a 
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+2L 2 {R p R a + P p P a ) - AQ 2 R p P a 
+gP° [-L 2 (R 2 + (FTL) 2 -Q 2 - L 2 )) 



(8) 



It is worth recalling that these expressions are obtained by anticipating the 
use of the relation 

L p P^ L (L) = (9) 

in order to drop any L p or L a in the expression of the Dirac's traces. Since 
this identity is not true for the gauge dependent part of the gluon propa- 
gator, one should not use these expressions of the traces to check the inde- 
pendence of the rate with respect to the gauge parameter £. Moreover, we 
discarded terms that will be killed later by the delta functions such as the 
one contained in S R (P) - S A (P) = 2ire(p )5(P 2 - M^). Since a 4-vector 
like P = (p , \Jj} 1 + M 2 ^ p) is not a linear function of the momentum P, we 
used some approximations to simplify the calculations, the effect of which 
is to neglect only terms that are always sub dominant 4 . 

We also notice that the statistical weights and delta functions present in 
Eqs. (3) and (6) are invariant under the change of variables P — > —R — L, 
L — ► L. Therefore, in the remaining factors of the integrand, we are allowed 
to drop the parts which are antisymmetric under this transformation. Col- 
lecting contributions from the two topologies, this symmetrization gives: 



Trace'" 7 1 . Traced 

| vertex , | self 

R 2 (FTL) 2 r 2 r 2 

n „ 2 (W P p 
-2Q 2 =- + 



sym 



2D 



R p R a + P p P° 



R 2 (P + L) 2 



R 2 (p + L) 2 I \R 2 (P + L 




pa 



-Q pU 1 



L 2 (L 2 + 



Q 2 L 2 



R 2 (P + L) 2 



+ 



+ (L 2 + Q ■ L) 



+ 



R2 (p + L) 2 



2 \(R 2 ) 2 ((P + L) 2 ) 2 J / 
which will be our starting point in the following sections. 

( Ml 



(11) 



4 For instance: 



R-P = Q(l + 



\p(p + q) 



(10) 
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3 Contribution to soft static photons 



3.1 Kinematics 

When the 3-momentum q of the emitted photon is zero, a lot of simplifi- 
cations occur. First of all, since there is one 3-vector less in the problem, 
we need only one angular variable which simplifies considerably the angular 
integrations. Moreover, as shown in [19], the non vanishing invariant mass 
of the emitted photon regulates all the potential collinear divergences when 
Q 2 /Qo ~ 1- Therefore, one can simply forget about the quark asymptotic 
thermal mass since the purpose of such a mass is precisely to regularize 
collinear singularities. This means that we can everywhere identify P and 
P at this level of approximation, since furthermore P is hard. 

From the identity S (P) — S (P) = 2ire(p )5(P 2 ), we extract the values 
p = ±p and r Q = q Q ± p. For the second cut quark propagator, we have the 
identity S H {R + L) - S A {R + L) = 2ne{r + l )5((R + L) 2 ) from which we 
extract the cosine of the angle 6' between p and /: 



Of course, we must require that this value be in [— 1, +1], which will re- 
duce the available phase space. This requirement leads to the following two 
inequalities: 



which lead to a phase space reduction that can be seen in Fig. 3, where 
the region excluded by the requirement cos 9' £ [—1,-1-1] has been shaded 
in dark gray. Other regions are excluded also by our choice of looking only 
at bremsstrahlung, i.e. excluding areas where L 2 > [25]. Finally, the only 
regions we have to consider are the unshaded ones. 

Having taken these constraints into account, the independent variables 
we are left with are for instance (the choice is not unique) r = p, l Q and I, ev- 
erything else being a function of these three. In particular, the denominators 
appearing in the calculations are 



cos 6' 



{r + l ) 2 -p 2 -l 2 
2pl 



(12) 



(l Q - I + p + q Q - p)(l + I + p + q Q + p) < 
(l Q - I + p + q Q + p)(l + I + p + q Q - p) > , 



(13-a) 
(13-b) 



R 2 = q {q ± 2p) 

(P + L) 2 = -q (q + 2l Q ± 2p) , 



(14-a) 
(14-b) 



the + sign corresponding to p Q = +p and the — sign to p D = —p. 
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Figure 3: Allowed domains in the (l a , I) plane for p = ±p. The area shaded in dark gray 
is excluded by the delta functions. The region shaded in light gray is above the light-cone 
(dotted lines). The solid curves are the transverse and longitudinal dispersion curves of the 
thermalized gluon. The vertical dotted line is the separation between e(p )e(r + 1 ) = +1 
and e(p )e(r + l ) = — 1. The value of the thermal mass has been exaggerated in order 
to make the figure more readable. 



It is worth examining more closely to which physical processes the regions 
I, II and III correspond. This is done just by looking at the signs of p Q and 
r + l . Examples are shown in Fig. 4. Besides the bremsstrahlung present in 
regions I and III, we see a new process in region II. This process corresponds 
to an annihilation of a quark- antiquark pair, one of the particles having 
previously undergone a scattering in the medium. Since the scale for the 
quark momentum is given by the temperature (p = 0(T)), and since we are 
looking here at soft photons, it is obvious from Fig. 3 that processes I and III 
have support of order T 2 in the (l Q , I) plane while process II has only support 
of order q a T, the integrand being the same. Therefore, we expect and we 
have checked that process II is suppressed by a factor of order q /T <C 1 
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Figure 4: Physical processes included in the diagrams of Fig. 2, in the region L 2 < 0. 
Region I: p < 0, r + l a < 0: bremsstrahlung with an antiquark. Region II: p < 0, 
r + l > 0: qq annihilation with scattering. Region III: p > 0, r + l > 0: bremsstrahlung 
with a quark. The particle on which the quark is scattered can also be a gluon. 



compared to bremsstrahlung. As a consequence, bremsstrahlung appears 
to be the dominant contribution as far as the I? < domain is concerned. 
In the remaining part of this section, we limit our study to the regions I 
and III. We can obtain a further reduction of the phase space by noticing 
that regions I and III give the same contribution since they are equivalent 
by a change of variables (indeed, after the symmetrization in Eq. (11), the 
integrands are invariant under the change of variables P — ► —R — L and 
L — > L). Physically, this means that photons are produced equally by 
quarks and by antiquarks. Therefore, we just consider region III {i.e. p Q > 
and r Q + l > 0) and multiply the result by an extra factor 2. Hence, the 
contribution of bremsstrahlung is given by 



AR N C f f f 

Imn Ai! /(g ,o) w 1 —^e 2 g 2 J dp J Idl J dl Q 

x \n F {r ) - n F (p )] [n B (l ) + n F (r Q + Q] 
Trace'" 7 



x £ Pa(L)P^(L) 

a=T,L 



Trace'* 7 1 
R 2 (PTL) 2 + R 2 R 2 



.(15) 



sym 
J Po=P 
cos«' = 



3.2 Extraction of the dominant terms 
3.2.1 General considerations 

Let us now concentrate on the matrix element to be integrated over the 
phase space depicted above. In order to perform the contraction of the 
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Dirac's trace obtained in Eq. (11) with the longitudinal projectors, we use 
the following expressions [27, 28] 



i> (L) = g p(J - U p U a + 



[L p - l U p )(L a - l U a ) 



I 2 



P P AL) 



-P P AL)+9 P *-^ 



(16-a) 
(16-b) 



where U = (l,o) is the 4-velocity of the plasma in its rest frame. An 
important simplification is obtained in the case of static photons since we 
have 

Q P P p a = • (17) 

Taking into account this simplification, we obtain after some algebra: 



]T Pa {L)P a pa {L) 



a=T,L 



Traced I . Traced ., 

I vertex I I self 

R 2 (FTL) 2 R 2 R 2 



sym 



(p T (L)-p L (L)) 
(Q + L) 2 



V(cosV-l) ( l2 



1 lL ^" /,'-'!/' • /.;-' 



t Q 2 {Q-L)^ 

'R 2 {P + L) 2 



+2 R 2 (P + L) 2 (QVW + ^VW) 
(Q ■ L) 2 



-2p T (L) 1-2 



R 2 (P + L) 2 



Q 2 L 2 



+ 



(R 2 ) 2 {{P + L) 2 ) 2 



(18) 



It is worth noticing that all the potential poles in l/L 2 (see Eq. (16-b)) 
have disappeared in this formula. This is nothing but a consequence of the 
gauge invariance of the set of diagrams we are looking at. 

Now, in order to extract the order of magnitude of each term in Eq. (18), 
we can use the following very rough rules: 

• R 2 ~ (P+L) 2 ~ 2pq Q which is correct for our purpose even if I is hard since 
I is bounded by a quantity proportional to p (see the previous paragraph on 
kinematics) . 

• n B (l ) ~ T/l Q which always gives the correct order of magnitude, even 
when l is hard. 

• n F (r + l ) = 0(l). 

• n F (r ) - n F (p ) w q n' F {p). 

• 1 - cos 2 #' = 0(1). 
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• The details of the x dependence are irrelevant to obtain the correct order 
of magnitude. 

• p TL (l,x) behaves like m 2 /(l 2 + m 2 ) 2 if one neglects its x dependence. 

AR 

Moreover, using the variable x = l Q /l, we can see Imll ^{q ,o) is a 
sum of terms like 

e 2 g\oJ dpn' F (p)pP j dlP J dxF(x)p T Jl,x) (y + 0(1)) , (19) 

where F(x) is a dimensionless function and where a + /3 + 7 = 3in order 
to give the correct overall dimension. A close inspection of Eq. (18) shows 
that a ranges from a = —1 to a = 3, taking all the integer values between 
these bounds. Eq. (18) shows also that (3 < 0. From the previous structure 
and since we can factorize a factor m 2 out of the spectral functions p T L , it 
is obvious that we can write the result as 



r AR. 



Imir'/(? ,o) « e 2 g 2 



m 2 T 



. P , * * 

\ m„ 1 



(20) 



where P(-) is a polynomial of two variables, of total degree 4. Under the 
assumption that q Q <C T, we are allowed to truncate this polynomial and 
keep only terms of degree in the variable q /T, leaving a polynomial of 
Qo/ m g only. Moreover, the above rules show that the term of Eq. (18) 
proportional to p 2 (cos 2 6' — 1)L 2 contributes to the constant term of this 
polynomial, and gives an integral that behaves like / dl/l for hard I. This 
means that a logarithm of order ln(l/<?) shows up in this coefficient. The 
argument of this logarithm can be written as T 2 /f(m 2 ,q 2 ) where /(•) is 
a function of dimension two. This function depends on both q a and m g 
since there can be a competition between q Q which appears as a kinematical 
infrared cut-off in the integral over dl and m g which appears in p T L and 
can also play the role of an infrared cut-off for the same integral. Using the 
same tools, it is quite easy to check that all the other coefficients of this 
polynomial are of order 1 (i.e. do not contain any large logarithm). 

Therefore, under the assumption that q a , m g <C T, we can formally put 
the result into the compact form: 



r AR 



ImU ^{q ,o)^Ce 2 g 



m 2 T 



2 2 «.* g 



q a 



\f{mi,q 2 ) \m gj 



(21) 



where C is a numerical constant and Q( ) is a polynomial. 
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3.2.2 Extraction of the logarithmic behavior 

Under some more restrictive assumptions, we can go further analytically. 
More precisely, it is possible to extract analytically the constant C in front 
of the logarithm of Eq. (21), as well as the function /(mg,g„) in the limit 
q <C m g . The assumption qo,m g <C T (i.e. g <C 1 and q Q <C T) ensures 
that the argument of the logarithm is large, so that the logarithmic term 
should be a fairly good approximation of the whole expression 5 . 

As mentioned before, the logarithm we are looking at comes from the 
term in p 2 (cos 2 9' — 1)L 2 in Eq. (18). Using Eq. (12), we obtain the following 
expression for the imaginary part of the photon polarization tensor: 

r, r, +00 +1 

imn /(g ,o)« J dp pip) J ~ { x ) 



o -1 



2p 



J dll 3 (p T (l,x)-p L (l,x)), (22) 



where 



-2Imn TL (x) 

Pt < l (/ ' X ' = (P(l - x 2 ) + Ren T L (x)) 2 + (Imn T L (x)) 2 (23) 

is the spectral function in the space-like region. This is the place where 
the condition q Q <C m g enters the picture. Indeed, if we don't make this 
assumption, the infrared regulator of the integral over dl will be a compli- 
cated combination of q Q contained in the lower bound and m g contained in 
the spectral functions. On the contrary, when q Q <C m g , only the largest 
regulator (i.e. m g ) plays a role, and the argument of the logarithm is quite 
simple. The integral over dl is elementary and yields an arctan function and 
a logarithm. Keeping only the latter 6 , we obtain: 



o 

x(Imnjx)-Imn T (x))ln(^j , (24) 



5 To summarize, the condition q ,m s <C T is essential to have a large logarithm, whereas 
the extra inequality q <C m g is necessary just to be able to calculate analytically the 
function f(m^,ql). 

6 The arctan term we discarded is convergent when performing the subsequent integra- 
tion over dx since the arctan is bounded by tt/2. 
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The terms neglected in that procedure show up only in the polynomial that 
would accompany the ln(l/g) in a more complete calculation, and are not 
tractable analytically. Moreover, for the same reason and because of the 
statistical weight in the dp integral that will cut off everything above p ~ T, 
we can replace the remaining logarithm 7 by ln(T A /m g ). Now, the dx and 
dp integrals are trivial and give: 

The production rate is then given by (see Eq. (2)): 

where the sum runs over the flavor of the quarks in the loop (e / is the electric 
charge of the quark of flavor /, in units of the electron electric charge). 

A comment is relevant concerning the sensitivity of the exchanged gluon 
to the hard scale. Indeed, the discontinuity of the effective propagator is 
used here in the space-like region, and the HTL approximation used to 
obtain this propagator may inaccurately reflect the phenomenon of Landau 
damping for a hard gluon. The consequence of this remark is that a loop 
correction on the gluon propagator may lead to an important three-loop 
correction to the photon emission-rate. 

Before comparing this analytical result with numerical estimates of the 
unapproximated expression, let us recall the domain in which this expression 
is expected to be valid. Firstly, we need the logarithm to be large in order 
to be dominant, which requires q ,m g <C T, i.e. j < 1 and q Q <C T. The 
additional purely technical condition is that q a <C m g , in order to keep 
simple the argument of the logarithm. On the plots of the figure 5, we 
show the ratio "numerical/theoretical", where "theoretical" is the formula 
given in Eq. (26) while "numerical" denotes a numerical evaluation of the 
contribution to bremsstrahlung of the complete matrix element as given in 
Eq. (18). The left plot shows that m g /T must be smaller than 0.1 in order to 
have an agreement between our approximations and the complete expression 
with an accuracy better than 5%. If m g /T is not small enough, then the 

7 It is possible to perform analytically the dp integral without this further approxi- 
mation. Doing so leads to a result in which the logarithm of Eq. (26) is replaced by 
\n(T 2 /m 2 g ) + 2(ln(7r/2) - 7) where 7 is the Euler constant. However, the additional 
constant is not complete since we have already neglected contributions to it in earlier 
approximations. 



dN 
dtdx 
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Fi gure 5: Comparison of numerical estimates of the complete matrix element with the 
simple theoretical expression obtained in Eq. (26). Both plots show the ratio "Numeri- 
cal/Theoretical". On the left plot, q /T is fixed at 1CP 4 and we look at the variations 
with rn g /T {i.e. with g). On the right plot, m g /T is fixed at 10 -2 and the photon energy 
varies from ultra-soft energies to hard ones. 



polynomial that comes with the ln(l/p) cannot be neglected anymore. On 
the second plot, we see that the approximations we performed inside the 
logarithm by assuming the smallness of q /m g are in fact valid far outside 
their expected domain of validity 8 , since we still have a reasonable accuracy 
with qo/m g ~ 10. 

3.3 Comparison with other approaches 

3.3.1 Extrapolation of the quasi-real soft photons results 

In a previous paper, we gave asymptotic formulae for the same quantity 
in the case where the photon invariant mass satisfies M^/T 2 <C Q 2 /q 2 < 
1 (see Eqs. (89), (90) and (93) of [19]). By extrapolating these estimates 
outside of their apparent range of validity towards the case of static photons 
for which Q 2 /q 2 = 1, we obtain exactly the formula of Eq. (26). Such an 
agreement means that, for a given energy q , the formulae established for 
the production of low invariant mass photons are very robust since they 
remain valid when extrapolated to the case of heavy photons at rest. 

From a technical point of view, this is made possible by the fact that the 
term that contains the collinear singularity in which we were interested for 
low mass photons and the term that develops the logarithm we extracted 
analytically for heavy photons are the same. 

8 This is presumably due to the fact that this extra assumption affects only the terms 
inside a slowly varying logarithm. 



17 



3.3.2 Braaten et al. results 



At the one loop level in the effective perturbative expansion, the produc- 
tion rate of soft static photons has been evaluated by Braaten, Pisarski and 
Yuan in [7]. The purpose of this paragraph is to present an analytic com- 
parison of their final result (Eq. (11) of [7]) with Eq. (26), in the domain 
q "C m g ~ gT <C T for which our expression has been justified. In this 
domain, we can retain from the result of [7] only the terms having the most 
singular behavior in q . Such terms are found only in the "cut-cut" part. 
Moreover, some of these terms develop a large logarithm ln(l/g) which is 
simple to extract analytically, in a way very similar to the method leading 
to Eq. (26). Applying these approximations to BPY's result leads to the 
following estimate for the one-loop production rate: 



dN 



dtdx 



_ p ^MX>?)(^) 4 >ng). m 



where vr? = is the soft quark thermal mass. Therefore, comparing 

with Eq. (26), we obtain the ratio: 



^Ibrcmss „ 32 N + N F /2 
^ll-loop ~ 3VT 2 C F 

which for 2 light flavors and 3 colors becomes 



(28) 



diV|t ™ s « ^ ~ 3.2 . (29) 

^ll-loop ^ 

This ratio is rather large, which means that bremsstrahlung is definitely 
an essential contribution to the soft static photon production rate by a hot 
plasma. 



3.3.3 Cleymans et al. results 

The bremsstrahlung production of a soft virtual photon has been considered 
in the context of the semi-classical approximation by Cleymans et al.. In 
their approach[12], they took into account the effect of the multiple scatter- 
ing of the quark in the plasma (Landau-Pomeranchuck-Migdal effect). In 
order to compare with our thermal field theory result we need to "undo" the 
effect of rescattering and consider only one collision of the (photon emitting) 
quark in the plasma. Cleymans et al. use several simplifying hypotheses: 
the energy of quarks or gluons is much larger than the temperature so that 
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Boltzmann distributions are used for particles entering the interaction re- 
gion and a factor 1 is assigned to those leaving it. The scattering of quark 
in the plasma is treated as in the vacuum, the only modification being the 
introduction of a phenomenological Debye mass m D to screen the forward 
singularity of the quark scattering amplitude. Neglecting furthermore the 
virtual photon momentum compared to the momenta of the constituents in 
the plasma the production rate can be factorized into a quark scattering 
term and a photon emission term so that the lepton pair rate can be written 
(see [19] for a similar expression in the case of a real photon) 



dN dq dq a r yr _ p . /T dp { dp 1 { 

dtdx ~ 24tt 4 QiJ .H 2 e (2^)32^^(2^)32^ 



x{2tt) a 5(P 1 + P 2 -P' 1 -P! 2 ) 

,\M\\P^P^)e> £ OU-^t m 

pol. e 

where \ftA\ 2 is the square of the matrix element of the quark scattering 
process. We have folded in the above expression the appropriate factor 
describing the decay of the virtual photon of mass Q into the lepton pair. 
With the above mentioned approximations and keeping the most singular 
term in the t— channel [12] we have 

\M\ 2 (P 1 ,P 2 ;P{,P^) = C ab 2g' (^^) (31) 

with C a b = 4/9 for quark-quark scattering and C ao = 1 for quark-gluon 
scattering. We then find for the rate of production of the pair at rest the 
following expression: 

where df = 2 S x 3 C = 6 and d = | 2x2jx2 s x3 c + 2 s x8 c = 26 + | are the 
degeneneracy factors introduced in [12]. Comparing with Eq. (26) we find 
for two light flavors 

| semi-class _ 15 



^-^Ibremss ^ 

i.e. the semi-classical result agrees with the thermal field theory result in its 
functional dependence but over-estimates the rate of production by about 
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50%. This difference appears to be due to the very approximate treat- 
ment of thermal effects of the dynamics of the plasma: for example, the 
ratio of quark-gluon scattering to quark-quark scattering is estimated to be 
N/ .5N F = 3 in thermal field theory (this is the ratio of the gluon contri- 
bution to the quark contribution in a hard thermal loop) compared to 1.5 
(taking account of the multiplicity factors associated to the quarks and the 
gluons) in the semi-classical approach. 

4 Contribution to hard real photons 

4.1 Kinematics 

Let us now concentrate on the case of hard real photons (T<q D = q). The 
kinematics for this situation is much more complicated than for static soft 
photons for two reasons. First, since Q is hard, we are no longer allowed 
to neglect Q in front of P or R. Moreover, since we are looking this time 
at real photons, we may encounter collinear divergences (like in the case of 
soft real photons [19]), and we must carefully keep the quark asymptotic 
thermal mass in the expressions. 

Now, from the identity S R (P)-S A (P) = 2ire(p )5(P 2 — M^) , we extract 
the value p a = ±o> p , where we denote u p = \J{p 2 + M^) and r a = q a ± ui p . 
The second delta function constraint provides us with the angle 9' between 
the 3- vectors r and I, via the relation: 

COBl) > = R2 - M ~y- r °'° +L2 . (34) 

Again, we must enforce the requirement — 1 < cos#' < 1, which implies the 
following set of inequalities: 

(r -r + l + l)(r + r + l -l) > (35-a) 
{ro-r + l - l)(r + r + l + I) < , (35-b) 

and leads to a reduction of the allowed domain in the (l ,l) plane. The 
above two inequalities may be rewritten as 

yj(r - I) 2 + Ml<l + r or l + r < - J ( r - I)* + M*, (36-a) 

-\A r + 2 + Ml<l + r < ^{r + I) 2 + Ml (36-b) 

The result of these inequalities is shown in Fig. 6, where the regions excluded 
by the requirement cos^' G [— 1, +1] are shaded in dark gray. In order 
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Figure 6: Allowed domains in the (l a , I) plane for p — ±lj p . The area shaded in dark 
gray is excluded by the delta function constraints. The areas shaded in light gray are 
above the light-cone (dotted lines). The light curves are the transverse and longitudinal 
mass shells of the thermalized gluon. The vertical dotted line is the separation between 
e(p )e(r + l a ) = +1 and e{p )e(r + l a ) = -1. 

to make the comparison easier with the case of static photons, we have 
also reproduced the boundaries of the allowed region for that case (i.e. the 
frontiers that one obtains in the limit M m — > 0). One can see that the old 
boundaries are asymptotes for the news ones. Again, we have three regions 
allowed by the above two inequalities for L 2 < 0. Again, regions I and 
III give the same contribution, as can be seen by the change of variables 
P — ► —R — L performed at an earlier stage of the calculation. From now on, 
we will drop region I and multiply by a factor 2 the contribution of region 
III. We start the discussion with the study of region III (bremsstrahlung 
processes) and then turn to region II (qq annihilation with scattering in the 
plasma, see Fig. 4). As will be seen, this region can no longer be neglected, 
contrary to the case of soft photon production. 
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4.2 Bremsstrahlung 

Since the delta function 5(P 2 — M^J makes more convenient the choice of 
p = \ \p\\ as an independent variable (instead of r), the quantities — r Q ± uj r 
that appear on the previous figure are to be understood as functions of p 
and the angle 8 between p and q. Therefore, this restriction of the allowed 
phase space is in fact a constraint that relates the independent variables 
p,0,l and l Q . In the following, doing the integration over 8 first will prove to 
be convenient since, due to collinear divergences, the result of this angular 
integral controls the order of magnitude of the result. If u = 1 — cos0 is 
the variable over which we integrate first, then the bounds on u depend on 
the other variables l,l Q and p because of the identities of Eqs. (36-a) and 
(36-b). In particular, the inequality Eq. (36-a) can be written in terms of 
u = 1 — cos 8: 

u < (p + q?-W{r + l o y-M^-lf ^ 
2pq 

which gives if one assumes that l and I are soft 

u<(^+ l ~yi-x), (38) 

where we denote x = l Q /l. As a consequence, since the relevant values of 
the momentum I are controlled by thermal masses of order gT, this upper 
bound is of order g. 

It is now worth giving expressions for the denominators that enter in 
the rate, since they are potentially dangerous when the photon is emitted 
collinearly to the quark. To that purpose, we need also another angular 
variable, which is left unconstrained by the previous considerations. This 
variable can be the azimuthal angle 4> between q and I when projected on a 
plane orthogonal to r. Therefore, if we denote 8" the angle between q and 
I, we can calculate this angle by 

cos 8" = cos 8 r cos 8' + sin 8 r sin 8' cos 4> , (39) 

where 8 r is the angle between r and q. This last angle is not independent 
of 8, to which it can be related by 

r sin 8 r = p sin 8 (40-a) 
r cos 8 r = p cos 8 + q ; (40-b) 
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so that Eq. (39) can be rewritten as 



cos 6" = - (cos cos 6' + sin 6 sin 0' cos <f>) + - cos 6' . 



(41) 



Using only the variables u = 1 — cos 6, p, I and x = l a /l, we can write [19] 



i? 2 - « 2pg 



M 2 
2p 2 



(42) 



and 



2?r 

/ 



dcp 



2vr(p + g) 



(P + L)2-M^ 



2qp 2 



1/2 ' 



(43) 



where we performed the integration over (j) at this stage since this denom- 
inator is the only place where (f> appears at the dominant order. As one 
can see, the first formula remains exactly the same as in the case of soft 
real photons [19], while the second one is only slightly modified by an extra 
(p + q)/p factor. As a consequence, the discussion made in [19] concerning 
the enhancement that one gets when performing the u integral is still valid. 
In particular, the terms where both denominators are present are enhanced 
by a factor of order p 2 /M^, while those where only one of them appear will 
get only a logarithm of this quantity. 

Again, our starting point is Eq. (11). The order of magnitude of each 
term is evaluated by taking into account the fact that the momentum Q is 
now hard, as well as the quark momentum. Moreover, we must take into 
account the enhancement by a factor of order 1/g 2 for terms having two 
denominators. Since the emitted photon is assumed to be real, we have 
Q 2 = 0. Therefore, one can check that only one term dominates in this 
matrix element: 



RPR 7 + PPP a 



R 2 (P + L) 2 ' 

which is again the same as in the case of soft real photons. 

Since the relevant values of u are of order M^/p 2 ~ g 2 , we have 



(44) 



r„ - iO r 



pq 

p + q 



u + 



M 2 

oo 

2p 2 



for p = +uj p 



(45) 



Therefore, one can see in Fig. 6 that the extra requirement L 2 < does 
not restrict significantly region III due to Eq. (45) (this is a result already 
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obtained in [19] according to which the kinematics confines the phase space 
to the L 2 < region in the collinear limit if L is soft). 

Since the term to consider is the same, and since the denominators have 
very similar expressions, it is completely straightforward to reproduce the 
calculation 9 that has been performed in [19]. The only minor difference lies 
in the contraction of R p R a + P p P a with the projectors, which gives a factor 
p 2 + r 2 approximated by p 2 + (p + q Q ) 2 in the collinear limit. Introducing 
the dimensionless quantities: 



w = w- (46 " a) 

~ ImII Ti (x) ~ ReII Ti (x) 

^0*0 = Tf 2 ' R tA x > = TH ' ( 46 " b ) 

lvl oo lvl oo 



we find 



T tt- 4 * uf ^ e 2 g 2 NC F . T T , T 
lmU /(Q) ~ y — A F (J T -J L )^ 



X 





i-oo 

J dp (p 2 + (p + q f) [n F (p) -n F (p + q )\ , (47) 



where we denote 
i 



j_ -. 

f dx ~ . , /" , Vw/u; + 4tanh + 4 

.]„,= ] — i„ . (x) dw — 

/ X T ' iV i (w + J R TL(a;) )2 + (/TL(x)) 2 



(48) 



As one can see, the integral over x and w can be factorized out of the 
expression and is the same as in the case of soft real photons [19]. In this 
previous work, we noticed that the integral over the gluon momentum could 
potentially be singular when the gluon is transverse, but is regularized by the 
quark thermal mass instead of the gluon thermal mass as one might expect. 
The interpretation of this result is in fact rather simple. Indeed, a close 
look at kinematics shows that the delta functions S(P 2 — M^) and S((R + 
L) 2 — M^) become incompatible in the limit L — > if > 0. Therefore, 
the quark thermal mass prevents an infrared divergence by reducing to zero 
the region of phase space where this divergence can occur. The only major 
difference with [19] concerns the integral over p, because of the fact that q Q 

9 As seen earlier, the kinematics gives an upper bound of order g for the variable u. 
Since we are considering enhanced angular integrals for which the relevant values of u are 
of order M^/p 2 ~ g 2 , this upper bound does not appear in the result of the u integral. 
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is now hard. Nevertheless, we can still give a rather compact expression for 
this integral in terms of poly-logarithms: 



dp(p 2 + (p + q ) 2 )[n F {p) -n F (p + q )\ 



3C(3) + y| + 

-\lo\/T, 



|) 2 ln(2) 



+4Li 3 (-e-^'/ J ) + 2| Li 2 (-e-l^l/ T ) 



-(f) ln(l + e-l^) 



where the poly-logarithm functions are defined via 

+00 n 
Li a(*) - £ -a ■ 



n=l 



(49) 



(50) 



Eq. (49) simplifies in the limit of extremely hard photons q S> T: 

+00 

/ dp{p 2 + {p + q ) 2 )[n F {p)-n F {p + q )] « q 2 T \n(2) . (51) 


Therefore, in this asymptotic regime, we find for the production rate of hard 
real photons: 

2^NC F aa s ( £ ej) ^~^ T (J T J L ) ln(2) , (52) 



dN 



dtdx 



bremss 9o2>T 



where e / is the electric charge of the quark of flavor / expressed in units of 
the electron charge. 

Besides this simple asymptotic result, it is worth adding that Eqs. (47) 
and (49) provide a generalization of the analogous formula of [19]. Indeed, 
the result we provide in the present paper for the bremsstrahlung contri- 
bution to real photon production is valid over the whole range of photon 
energies, and in particular reduces to Eq. (55) of [19] in the limit of soft q Q . 

4.3 qq annihilation with scattering 

The discussion of the contribution of region II can be carried out in a similar 
way. We recall that p Q = —lo p < so that Eq. (42) becomes 



R 2 -M 2 



-2pq 



v + 



M 2 
2p 2 



(53) 
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with the notation v = 1 + cos 9. For hard enough q Q , the statistical weight 
n F (Po) ~ n F (Po + Qo) is equal to 1 for — q Q < p Q < and equals everywhere 
else 10 . Therefore, we can restrict to the range < p < q Q and we check that 
Eq. (43) remains valid except for the changes u — > v and q + p — > q — p. The 
enhancement mechanism in the terms that contain the two denominators 
goes through as before due to their behavior near v ~ g 2 . 

Such terms are dominated by the region where v ~ g 2 , which enables us 
to obtain — r Q — oj r w —2{q —p) and — r + oj r ~ 0{g 2 T). As a consequence, 
the boundaries in the (l , 1) plane are — 1 < x < 1 and < / < 2(q—p)/(l—x), 
where x = l a /l. Therefore, the integral over I and x gives the same factor 
J T — J L as before. The integral over p in Eq. (47) is now to be replaced by 

j 'dp(p 2 + (q -p) 2 ) = ^, (54) 


which leads to the following asymptotic contribution for the region II: 

^ 1% N C F aaJye 2 f )^e^°/ T (J T -J L ). (55) 
n <? »T3 7r 5 F s \y } > q Q Ll y ' 

Therefore, it appears that region II dominates over bremsstrahlung in the 
asymptotic regime. 



dN 
dtdx 



4.4 Comparison with previous results 

The production rate of hard real photons has already been calculated at the 
one-loop order in [10, 11] as an application of the effective theory based on 
the resummation of hard thermal loops. We now compare the contribution 
of bremsstrahlung obtained above with their result. For hard real photons, 
the predictions of [10] are 



dN 



dtdx 



l-loop <?o>T 



0^W£4)^ ln (.^), m 



with the constant c ~ 0.23. 

It is worth recalling here that the quantities J T L that appear in Eqs. (52) 
and (55) are functions of the ratio M^jm^ when Q 2 = 0, i.e. depend only 

10 The transition between and 1 takes place in a range of width T in the variable p . 
These side-effects are negligeable if the condition q ^> T is satisfied since they provide 
corrections of relative order T/q . 
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on N and N F . For N = 3 colors and N F = 2 light flavors, we can eval- 
uate numerically J T ps 4.45 and J L —4.26. The following plot shows 
a comparison of bremsstrahlung, c/g annihilation with scattering, and one- 
loop contributions. In this figure, the bremsstrahlung contribution is taken 

• * fane-loop □ 

1.4 - o • Bremsstrahlung o 

• . q-q with scattering a 

o . . Total . 



1.2 - 
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log r»,/r; 
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Figure 7: Comparison of various contributions to Imll^(Q) for a hard real photon. The 
comparison is made for N = 3 colors and N F = 2 flavors. The value taken for the coupling 
constant is g — 0.5. 

from Eq. (47); for the region II, we use the asymptotic result obtained for 
large q Q 11 , and the one-loop result is a numerical evaluation of the dia- 
grams considered in [10]. The value used here for the coupling constant is 
g = 0.5. We can see that the bremsstrahlung is the dominant contribu- 
tion for the smaller values q a , whereas the region II becomes dominant for 
hard enough q Q . For intermediate photon energies (around T), the three 
contributions have equivalent orders of magnitude. For higher values of the 
coupling constant g, the relative importance of the one-loop contribution 
tends to decrease. We can also see that the sum of the three contributions 
is significantly above the one-loop contribution considered alone. 

n This asymptotic estimate may be incorrect for smaller values of q . Nevertheless, we 
expect it to decrease like q since its support in the (l ,l) plane decreases like q , so that 
region II is certainly subdominant for q a < T. Therefore, the asymptotic result is sufficient 
in this comparison. 
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5 Conclusions and perspectives 



We have discussed the production of a real photon or of a lepton pair in 
a hot quark-gluon plasma in thermal equilibrium. We assumed that the 
plasma is described by thermal Quantum Chromodynamics and we worked 
in the framework of the effective theory obtained after resummation of hard 
thermal loops. We have shown that, at leading order in gT, it is necessary 
to include the one-loop as well as the two- loop diagrams. The one- loop 
diagrams of the effective theory correctly account for the contribution of 
soft fermion momenta in the loop, whereas in the two- loop diagrams hard 
momenta play a dominant role. 

Many physical processes are contained in the effective theory calculated 
up to two-loop order. To simplify, it can be said that, at one-loop, scattering 
processes mediated by (soft) fermion exchanges are of paramount importance 
for the emission of the photon. The corresponding photon production rate 
consequently strongly depends on the soft fermion thermal mass m F . On 
the other hand, the two-loop topologies account, among other possibilities, 
for bremsstrahlung processes where the photon is emitted by a hard quark 
scattering in the plasma via a (space-like) gluon exchange: in that case the 
rate is proportional to the square of the gluon mass m g . We have calculated 
the contribution of such processes. 

For the production a soft virtual photon (q Q ~ gT, q = o) the bremss- 
trahlung processes largely dominate over the one-loop result in the range 
Qo < m F • For the case of real photon production an interesting enhancement 
phenomenon occurs in the bremsstrahlung processes. Because of the van- 
ishing photon mass, the fermion propagators become infinite when a quasi 
forward-scattered quark emits a collinear photon: the singularities are reg- 
ularized by an interplay between the thermal masses of the fermion and the 
gluon. It was seen before that this leads to an enhancement by a factor 1/g 2 
of the bremsstrahlung contributions to the real soft photon rate so that the 
two-loop diagrams entirely dominate over the one-loop contribution. On the 
other hand, the bremsstrahlung production of a hard real photon occurs at 
the same order in g as the one- loop result. A rather simple analytic ex- 
pression, valid for soft, hard and ultra-hard photon has been derived. For 
ultra-hard photons, another process becomes dominant, consisting of qq an- 
nihilation where one of the fermions undergoes a scattering in the medium. 
In all cases, the calculated two-loop contributions considerably increase the 
rate of photon or lepton pair emission. Our results can easily be numeri- 
cally extended to cover, on the one hand, the case of a soft lepton pair at 
non-vanishing momentum, and on the other hand, the case of lepton pairs 



28 



produced at large momentum with a small invariant mass. All the features 
of the present results should survive in these more general situations. 

A word of caution should be given concerning the bremsstrahlung con- 
tribution to the soft virtual photon rate. The result shows a sensitivity to a 
hard space-like gluon and therefore one may suspect that the extrapolation 
of the effective gluon propagator in the hard region is not complete: to be 
consistent may require taking into account three-loop diagrams in very much 
the same way two-loop diagrams were needed besides the one-loop diagram 
with a hard space- like fermion propagator. On the contrary, for the case 
of a real photon emission, the bremsstrahlung rate is sensitive only to soft 
gluon momenta, and the calculation is therefore expected to be complete. 

Our study does not cover all the physical processes included in the two- 
loop diagrams. In particular, the contribution with a time-like gluon should 
be added. However such processes when the gluon is hard are already con- 
tained in the one-loop diagram: counterterms should therefore be included. 
To do this will be a very interesting practical exercise in the use of the 
effective theory up to two- loop. 
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A Hard thermal loops and counterterms 

When using effective theories based on the summation of hard thermal loops 
at higher orders, there is potentially a possibility to have multiple counting 
of thermal corrections that should be there only once. 

As a first illustration of this problem in a trivial context, let us use 
the example of a massless real scalar field with a A 2 (/> 4 interaction in 4 
dimensions. The Lagrangian of such a model is 

£=±0^0-^. (57) 

After the calculation of the one-loop tadpole, one realizes that this diagram 
- a hard thermal loop in the terminology of [3] - generates a thermal mass 
m = AT-y/24 that can be important for the phenomenology of soft modes. 
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Therefore, the idea of the HTL resummation is to include this thermal mass 
in an effective Lagrangian 

£ eff = \d»W* ~ \m 2 4> 2 - |V , (58) 

which is the Lagrangian of a real scalar field of mass m. 

Let us assume now that one uses this effective theory instead of the bare 
one to calculate the same tadpole diagram. The result would be 

„ A 2 T 2 / 3m 3 /m\ 2 fm\ \ 9 . . . 

11 " ~W (» — T ~ 2? ( t) ( t) + ■ ■ ■ J " ™ (1 + ° (A)> ■ <59) 

and since we start now from a propagator with squared mass m 2 , the resum- 
mation of the self-energy II would lead us to a propagator with a squared 
mass 2m 2 (1 + (9(A)), approximately twice larger than the correct thermal 
mass. Obviously, the above result arises due to multiple countings of the 
same thermal correction. Stated differently, this is a consequence of the fact 
that this effective theory is more than a mere reordering of the perturbative 
expansion of the bare theory since its Lagrangian is different. To solve this 
problem, one must write 

C = C cS + Cct with = ^m 2 (f) 2 , (60) 

and treat the counterterm £ ct as an interaction term, just like A 2 </> 4 /4!. The 
effect of this counterterm is of course to subtract at higher order the thermal 
corrections that had already been included at the tree level via the effective 
Lagrangian, in order to avoid multiple countings. For instance, in the above 
example, the correct answer for the tadpole II when one takes care of the 
counterterm is 

„ A 2 T 2 / 3m 3 fm\ 2 ^ I ' m\ \ 2 „^ 2 , . . 

n " ~u i 1 - H t - (r) (r) + " j - m = 0[Xm ) (al) 

which is a perturbative correction to the thermal mass found at the previous 
step, as it should be. When this counterterm is taken properly into account, 
the effective theory is nothing more than a reordering of the perturbative 
expansion, since the overall Lagrangian remains unmodified. 

This problem also arises in effective gauge theories, where the situation 
is a bit more complicated since we need there an infinite series - one for 
each hard thermal loop - of non-local counterterms. These counterterms 
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are defined to be the opposite of the HTL contribution to the corresponding 
function. Then, to a given diagram obtained in the effective theory, one 
should add the diagrams obtained by collapsing each loop in turn and re- 
placing it by the corresponding counterterm. In order to be more definite, 
let us consider the example of the two- loop diagram of figure l-[b]. In the 
figure 8, we have represented next to the one-loop contribution to the po- 
larization tensor of the photon some of its two-loop corrections. It should 




Figure 8: Some corrections to the one- loop contribution (left diagram) to the photon 
polarization tensor. A black dot denotes an effective propagator or vertex. Crosses are 
HTL counterterms. The boldface loop may reproduce what is already included in the 
one-loop diagram when its momentum is hard. 

be clear from the figure that when the boldface loop of the second diagram 
carries a hard momentum, then this loop reproduces the HTL part of the 
effective vertex in the dotted box already included at one-loop, plus new 
sub-dominant perturbative corrections. The purpose of the third diagram is 
precisely to subtract a quantity equal to the HTL contribution to this ver- 
tex, so that what remains constitutes only new contributions. The net effect 
of this procedure is thus to reorder the terms of the perturbative expansion. 

Let us now explain why, despite their conceptual importance, the coun- 
terterms do not contribute in the case of bremsstrahlung. As said in section 
2, the contributions of bremsstrahlung to photon production come from 
two- loop diagrams in which the gluon propagator is cut, and where one 
retains only the Landau damping part (L 2 < 0) of the cut gluon propagator. 
Technically, the gluon included in the 3-point HTL that appears on the third 
diagram of Figure 8 is a bare one, and therefore its discontinuity has support 
totally included in the time-like region. As a consequence, the diagram with 
counterterms contributes only in the region of phase-space where I? > 0, 
and cannot contribute to bremsstrahlung. Even if they are not a worry 
in the case of bremsstrahlung, it was important to discuss the potential 
effect of counterterms since, having shown that the second diagram gives an 
important contribution, it is not enough to conclude that the two-loop order 
is important, as it may be canceled by counterterms. 
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B Phase space considerations 



The purpose of this appendix is to emphasize the importance of being in a 
thermal bath in order to have a hard phase space for the quark loop. This 
is indeed crucial in the calculations performed in the previous sections since 
this feature of the diagrams considered in this paper enables them to be of 
the same order of magnitude as one-loop diagrams. In order to make the 
following discussion more intuitive, let us first transform the bremsstrahlung 
photon production rate in a way that separates more clearly the phase space 
from the amplitude of the process producing the photon. The tools to do 
that have already been presented in [19] (see the section devoted to the 
comparison with semi-classical methods), so that we only give the result 
here 12 : 



where R = P + Q. This formula tells us that in order to calculate the 
contribution of some process to the photon production rate, we just have 
to integrate the amplitude squared of this process over the phase space of 
unobserved particles (here, the incoming and outgoing quarks). When doing 
this integration, the external particles are put on their mass shells, and 
accompanied by the appropriate statistical weight. The main advantage 
of this formula is that it exhibits a clear separation in two factors: the 
amplitude squared of the process producing the photon, and the phase space 
of the quarks. Therefore, this formula can be used to separate the estimation 
of the order of magnitude of the diagram in two steps: the order of magnitude 
of the amplitude, and the size of the phase space over which it must be 
integrated. 

Looking at this formula, it is clear that the effect of the thermal bath 
appears only in the phase space. Indeed, the amplitude that appears in 

12 This formula is given here for the production of real photons by bremsstrahlung. For 
other processes, analogue formulae still hold, in which the amplitude squared have to be 
appropriately modified. 




X t + 



— > 1 > K+L > * > 



x2tt5(P 2 - Ml) 2nS((R + L) 2 - Ml) 

x2tt5(K 2 - Ml) 2tt6((K + L) 2 - Ml) 

xn F (r D + l )n F (k )[l - n F (p )][l - n F (k + l Q )\ , 



(62) 
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Eq. (62) is nothing but a zero temperature one (it does not contain any 
statistical weight). Stated differently, if, instead of looking at photon pro- 
duction by a thermal bath, we were looking at bremsstrahlung photon pro- 
duction in p-p collisions (i.e. the two scattering quarks come from protons 
in colliding beams) , the statistical factors of incoming quarks would have to 
be replaced 13 by structure functions 14 of a quark inside a proton, and the 
statistical weights of outgoing quarks would be replaced by 1. This difference 
is precisely the point which makes the thermal bath dramatically different 
from the p-p collision. Indeed, in the case of the thermal bath, the statisti- 
cal functions have a support which is 3-dimensional since the plasma is an 
isotropic medium. On the other hand, the structure function of the quark 
inside protons of the beam is vanishing if the quark has a direction different 
from that of the beam. This will make a difference when one performs the 
integral over the momentum p. Indeed, if one has something like / dpp n f(p) 
in the case of the thermal bath, the integral would be / dpp n ~ 2 f(p) for the 
case of p-p collisions. As a consequence, the integral over the quark mo- 
mentum is more likely to be sensitive to hard momenta in the case of the 
thermal bath. Therefore, even if it is not a rigorous proof, these consid- 
erations show why some higher order processes which are not dominant at 
zero temperature may become dominant in a plasma, due to a bigger size 
for their phase space. 

Moreover, Eq. (62) may help to understand why two-loop contributions 
may be as important as one-loop ones. Indeed, it shows that the final or- 
der of magnitude of a contribution results from a competition between two 
effects. The first one is the order of magnitude of the amplitude, which 
usually becomes smaller when the number of loops increases since the num- 
ber of coupling constants increases also. The second aspect of the problem 
is the size of the phase space, since it may happen that due to kinematical 
constraints, the one-loop phase space is much smaller than the two-loop one. 
Both effects can compensate so that two-loop diagrams contribute also at 
the dominant level. This is precisely what happens in the case of photon 
production by a plasma: the one-loop phase space is soft due to kinemat- 
ical constraints, while the two-loop thermal phase space can be hard (as 
explained above, this is possible because we are in an isotropic medium). 

13 To understand the analogy between the two situations, one may see a proton beam as 
a dense medium containing quarks and gluons, with distributions related to the structure 
functions of quarks and gluons inside the proton. Such a medium is highly anisotropic, 
since all the partons go in the direction of the beam. 

14 More exactly, to obtain the correctly normalized photon production rate, one should 
also multiply by the proton densities in the beams. 
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Therefore, the smallness of the two-loop amplitude is compensated by the 
size of the two-loop phase space. 
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